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Abstract. The wavelength dependences of interstellar extinction and polar- 
ization, supplemented by observed elemental abundances and the spectrum of 
infrared emission from dust heated by starlight, strongly constrain dust models. 
One dust model that appears to be consistent with observations is presented. To 
reproduce the observed extinction, the model consumes the bulk of interstellar 
Mg, Si, and Fe (in amorphous silicates), and a substantial fraction of C (in car- 
bonaceous material), with size distributions and alignment adjusted to match 
observations. The composition, structure, and size distribution of interstellar 
grains is the result of injection of dust from stellar outflows into the interstellar 
medium (ISM), followed by destruction, growth, coagulation, and photoprocess- 
ing of interstellar grains. The balance among these poorly-understood processes 
is responsible for the mix of solid material present in the ISM. Most interstellar 
grain material present in the diffuse ISM must be grown in the ISM. The amor- 
phous silicate and carbonaceous materials that form the bulk of interstellar dust 
must therefore be the result of grain growth in the presence of ultraviolet radi- 
ation. Dust in high-z systems such as J1148-I-5251 is also produced primarily 
in the ISM, with supernova-produced dust contributing only a small fraction of 
the total dust mass. 



Introduction 



Trumpled (|l930l ) and iLindblaJ (|l935l 1 recognized that reddening of distant stars 



was produced by submicron particles present in interstellar space, but over 70 
years later we continue to work to determine the composition, shape, and size 
distribution of interstellar dust. Despite substantial advances in our under- 
standing of interstellar physics, the problem is far too difficult for an ah initio 
approach. We do not know enough about physical conditions in stellar out- 
flows (where some grains are formed) and the interstellar medium (where grains 
undergo modification) , and our understanding of the physical and chemical pro- 
cesses occuring on grain surfaces and within grains is very limited. We cannot 
predict with any confidence how grain materials will evolve over hundreds of 
Myr in the interstellar medium (ISM). We need to use what physics we know, 
but observations must guide us. 

Almost all of the observational information concerning interstellar dust 
arises from interaction of electromagnetic radiation with dust: scattering, ab- 
sorption, and emission . The fir s t "de tection" of interstellar dust was through 
extinction of starlight: iHerschel ( 17851 ) remarked on the absence of stars in cer- 



tain portions of the Milky Way ( "An Opening in the Heavens'''), although he did 
not understand its cause. Over a century later, iBarnar 3 (ligOTl ) noted that the 



neighborhood of p Oph was an example of "an apparently absorbing medium". 
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Figure 1. Average extinction vs. for dust in interstellar diffuse clouds. 
For > ljim~^, this is based on the parameterization of extinction by 
iFitzpatricS (Il999f). with addition o f the 7 7 strongest diffuse interstellar bands 
(DIBs) from iJenniskens fc Desert ()1994|) . Three of the strongest DIBs are 
labelled. 



We now study scattering, absorption, and emission of radiation by interstellar 
dust at wavelengths ranging from A (X-rays) to cm (microwaves). 



2. Observational Constraints 



2.1. Extinction 

The wavelength-dependent extinction of starlight - the so-called "extinction 
curve" - remains the principal source of information about interstellar dust. 
Interstellar extinction has been studied on many sightlines at wavelengths be- 
tween about 2//m and 0.1/im. On sightlines with sufficient column density, it is 
possible to study the extinction at wavelengths as long as ~ 20/im. An average 
extinction curve is shown in Figure [TJ 

The extinction curve contains spectral features that constrain the compo- 
sition of the dust. The strongest feature by far is a broad "bump" peaking 
near 2175A. The strength of this feature requires that it be produced b y a sub- 
stanc e composed of high-abundance elements, such as C, Mg, Si, or Fe ( Draind 
19891 ). The position of the feature, and its width, are strongly suggestive of 
TT — > TT* excitations in aromatic car bon, such as g raphite or polycyclic aromatic 
hydrocarbons. Some authors (e.g., iDraine Sz Lil 1200 7) think that the feature is 
produced by the large population of polycyclic aromatic hydrocarbons that is 
required to explain a number of infrared emission features. 
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The extinction curve in Fig. [T] shows a conspicuous feature near 10/im that 
can be confidently attributed to Si-0 stretching modes in amorphous sihcates; a 
companion feature at 18//m is attributed to 0-Si-O bending modes. The absence 
o f fine structure in the profile indicates that the silicate material is amorphous 
- Kemper et al. ( 20051 ) place an upper limit of 2% on the crystalline fraction. 



There is a weak absorption feature at 3.4/xm that is attributed to the C-H 
stretch in hydrocarbon material. The nature of th e hydrocarbon material is 



controversial: from the profile shape and strength, iPendleton Sz Allamandola 
(|2002i ) estimate t hat the carbonaceou s material is ~85% aromatic and ~15% 
aliphatic, whereas Dartois et al.l ( 20041 ) conclude that at most 15% of the carbon 



is aromatic. 

Finally, we note that the extinction curve in Fig. [1] actually contains ~200 
weak but detected spectral features - the "di ffuse interste llar bands", or DIBs. 
86 years after the first DIBs were reported bv iHeger ( 19221 ) . it is embarassing to 



have to admit that not a single one has yet been identified. 

The extinction curve shows a smooth rise from the near IR to the vacuum 
UV. A broad range of grain sizes is required to reproduce this: the rapid rise in 
extinction shortward of 1500 A requires very large numbers of grains with radii 
a < 0.02 /im, but the extinction in the visible requires that most of the grain 
mass be in grains with sizes 0.05//m < a < 0.3/xm. The overall strength of 
the extinction per H nucleon requires that the grains providing the bulk of the 
extinction be based on high-abundance elements such as C, O, Mg, Si, and Fe. 

Extinction curves vary from one sightline to another. Figure [1] shows an av- 
erage extinction curve, but grain models must be able to account for the observed 
sightline-to-sightline variations in the wavelength-dependence of extinction. 

2.2. Other Evidence 

In addition to the wavelength-dependent extinction, a modern dust model must 
also be consistent with a wide range of other observations of electromagnetic 
radiation absorbed, scattered, or radiated by dust: 

• Polarization of starlight by interstellar grains. 

• Scattering of starlight by grains (reflection nebulae, diffuse galactic light). 

• Small-angle scattering of X-rays by interstellar grains. 

• Infrared and submm emission from interstellar grains heated by starlight. 

• Polarization of the infrared and submm emission from dust. 

• Microwave emission from dust. 

In addition to the abov e, we can study the i n terstellar grains tha t are entering 
the heliosphere today ( Landgraf et al.l [20001 : iKriiger et al.l [20071 ) . and we can 



analyze presolar grains trapped 4.5 Gyr ago in primitive meteorites and comets. 

Our estimates of the total interstellar abundances of elements such as C, Mg, 
Si, and Fe, plus direct measurement of gas-phase abundances using absorption 
lines, tell us how much of each element appears to be locked up in solid grains. 
Gas-phase abundances of species such as Si, Ca, and Ti are observed to vary 
from one sightline to another, which implies that depletion from the gas must be 
able to occur relatively rapidly; a successful grain model must provide enough 
surface area for depletion of Si~^, Ca"'", and Ti"'" to take place rapidly. 
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The observed sightline-to-sightline variations of gas-phase D/H is best un- 
derstood if D is actuahy de pleting onto dus t grains (Draine 2004, '2006a) and evi- 



dence for this is now strong (jProchaska et al . 2005; Linsky et _al. 2006: Lahement et al 



20081 ). A dust model must contain a component capable of trapping a substan 



tial f raction of the D present in the ISM. The PAHs may be able to accomplish 
this (|Draindl2006al ). 

In general, dust grains play a crucial role in interstellar chemistry, and a 
successful dust model should eventually be able to "predict" the rate for forma- 
tion of interstellar H2 via reactions on grain surfaces. Photoelectrons emitted 
by interstellar grains make a major contribution to heating of interstellar gas, 
and a successful grain model must be in quantitative agreement with heating 
rates inferred from observations of interstellar gas. 



3. A Model Using Amorphous Silicate and Carbonaceous Grains 

Current modeling techniques do not allow us to "invert" the observational con- 
straints to arrive at a grain model. Constrained by elemental abundances, all 
sensible grain models consist predominantly of carbonaceous and amorphous sil- 
icate material, but there are various candidates for the carbonaceous material. 



For example. IZubko et al. ( 20041 ) consider models involving varying amounts of 



PAHs, graphite, 3 different types (AC, BE, ACH2) of amorphous carbon, and 
an organic refractory material with composition C25H25O5N. 

Here we will focus on one particular model, in which the carbonaceous mate- 
rial in ultrasmall particles is assumed to have the optical and thermal properties 
of PAHs, while the larger carbonaceous grains are assigned the optical properties 
of graphite. If polarization is no t of interest, the particl e s may be approximated 
as spheres. For these materials, IWeingartner Sz Draine ( 200ll . hereafter WDOl) 



found size distributions that reproduce the extinction in various regions of the 
Milky Way, and in the Large and Small Magellanic Clouds. For current es- 
timates of solar and interstellar abundances, the WDOl size distributions use 
somewhat more C, Mg, and Si than is nominally available, but one should keep 
in mind that (a) there are modeling uncertainties (e.g., the grains are certainly 
not solid spheres), (b) the actual interstellar abundances of C, Mg, Si, and Fe are 
uncertain, (c) gas-phas e abundances of C have recently been revised downward 
(|Sofia Parvathill2009l ). increasing the amount of carbon inferred to be in dust: 



it now appears that only ~ 1/3 of the C is in the gas, and 2/3 is in the dust. 

The PAH-graphite-amorphous silicate model satisfactorily reproduces in- 
terstellar extinction (see WDOl), but it can also be tested in other ways. A 
dust grain illuminated by starlight will absorb energy, which will be reradiated 
in the infrared. The infrared emission spectrum will depend upon both the 
temperature of the grain and its infrared opacity. 

Large a ^ 0.02/im grains in the local starlight background will be heated 
to a more-or-less steady temperature of 15-20K. However, very small grains (1) 
absorb photons much less frequently, and (2) have very small heat capacities, so 
that one absorbed photon can raise the grain to a high temperature, followed by 
very rapid cooling. Because of this, one cannot speak of an average temperature 
for very small grains - instead, one must determine the temperature distribu- 
tion function (dP/dT), where dP is the probability of finding the grain in the 
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Figure 2. Global emission from NGC 7331 (from|Draiiie et al."2007'). Sym- 
bols are 2MASS, Spitzer, IRAS, and SCUBA photometry. The sohd curves is 
generated by the DL07 dust model, which requires a dust mass 1.1 x IO^Mq 
to reproduce the observed emission spectrum. 



temperature interval [T, T + dT] . The infrared emission is given by 

#=E T.?fdT(§,) 5.(r)a,.,,(A) . (1) 



compi size J - 

To calculate the emission from a grain model that includes very small grains, 
one must first find {dP/dT)ij for the different compositions i and grain sizes j 
that are present in the model. For large grains, dP/dT 6(T — (T)), where 6 
is the Dirac delta function, but for small grains dP/dT has a tail extending to 
high T. 

The DL07 grain model uses the WDOl size distribution for all except the 
PAHs; the size distribution of the PAHs was adjusted to try to reproduce the 
A < 25/im infrared emission. DL07 have solved for the temperature distribution 
functions dP/dT for different grain sizes and composition, and have calculated 
the resulting infrared emission spectrum for various intensities of the starlight 
heating the grains. 

The diffuse ISM away f rom the Galactic plane in the Milky Way has been 
measure d by COBE-DIRBE (iDwek eFaPliggTl - lArendt et aUfTooa ^ and by COBE- 



FIRAS (|Finkbeiner et alJ Il999l 1~ The emission predicted for the DL07 grain 



model is in good agreement with the observed diffuse emission from high-latitude 
regions in the Milky Way. 

If we assume that grains in other galaxies are like grains in the Milky Way, 
we can test the grain model by comparing it to the emission observed from other 
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Figures. Emission fro m centr al regions of 6 SINGS galaxies. Error bars are 
for 5-33^m IRS spectra (jSmith et al. 2007). Rectangles are IRAC and MIPS 
photometry. Solid line is a model consisting of starlight plus emission from 
dust heated by starlight. The DL07 physical grain model is able to reproduce 
both the PAH features and the IR continuum for diverse galaxies, but with 
different PAH fractions gpAH and different starlight intensity distributions. 
From .Drainc ct al.. (200a) . 
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galax ies. The Spitzer Infrared Nearby Galaxy Survey (SINGS, iKennicutt et al 



I2OO3I ) consists of 75 nearby galaxies observed with all of the inst ruments on 
Spitzer Space Telescope (Werner et al. 2004). iDraine et al.l ( 20071 ) have used 
the DL07 dust models to estimate the dust content and starlight intensities for 
65 of the SINGS galaxies. Figure [2] shows the observed SED for one example, the 
well-observed Sb galaxy NGC 7331. With suitable assumptions for the starlight 
heating the dust, the DL07 dust model closely reproduces the observed infrared 
and submm emission. 

Modeling the IR emission from NGC 7331 yields an estimated dust mass 
Mdust = 1-1 X 10^ A/0. The HI 21cm emission give an HI mass Af(HI) = 1.0 x 
IO^Mq, and the CO 1^0 luminosity gives an estimated H2 mass M(H2) = 
1.6 X 1O^°M0 (for an assumed CO "X factor" 4 x lO^^cm"^ K km s"^), resulting 
in a dust-to-H mass ratio Mdust/A^H = 0.0043, reasonably close to the expected 
value ~ 0.007, especially when one considers the uncertainties in estimation of 
the H2 mass from the CO 1^0 luminosity. 



4. Polarization of Starlight and Polarized IR-Submm Emission 

WDOl showed that a model based on amorphous silicate and carbonaceous 
grains, including PAHs, could reproduce observations of extinction in the Milky 
Way, LMC, and SMC, and DL07 showed that these same dust models could 
reproduce observations of infrared emission. 

In addition to interstellar extinction, a dust model should be able to re- 
produce observations of polarization of starlight by aligned dust grains. The 
infrared emission from aligned dust grains will be polarized, and microwave and 
submillimeter polarimetry has advanced to the point where it is beginning to be 
possible to measure the polarized emission from the diffuse ISM. For computa- 
tional convenience, the WDOl and DL07 models assumed spherical grains, and 
therefore couldn't say anything about polarization. 

Early modeling of starlight polar ization at optical wavele ngths approxi- 
mated the grains by infinite cylinders ( Greenberg Sz Hone 197 j ) , but improved 



numerical methods and faster computers allowed Kim &: Martini ( 19951 ) to model 



the grains as spheroids. Advances in polarimetric techniques at submm to cm 
wavelengths now allow detection of polarized emission from the diffuse ISM (e.g. 
iPage et al.ll2007l ) and theorists should now make predictions for the degree of 
polarization that is expected for different dust models. 

In addition to reproducing the wavelength-dependent extinction, a dust 
model must first reproduce the polarization of starlight as a function of wave- 
length. All dust models must have aligned, nonspherical grains containing amor- 
phous silicates, because the interstellar sil icate feature is known to produce po- 
larization in extinction ( Smith et al. 2000l ) . If carbonaceous material is present 



in a separate grain population, at this time there is no evidence that such grains 
are aligned. Absence of polarization in the 3.4;um C-H stretch feature suggests 
that the carbonaceo us particles (at lea. st those responsible for the 3.4/im feature) 
may not be aligned ( Chiar et al.ll2005l ). 



Using models based on carbonaceous grains (including PAHs) and amor- 
phous silicate grains with various assumptions about (1) the sh apes of the grains 
and (2) whether or not the carbonaceous grains are aligned, IDraine &: Fraissi 



8 



B. T. Draine 




Figure 4. Left: Infrared emission vs. wavelength A for 4 models (modl- 
mod4) with spheroidal grains that reproduce extinction and polarization 
of starlight (jPraine fc Fraiss3 l2009( l and for the spherical grain model of 
iDraine &: Lil |2007, DL07). Crosses show the ernission per H for loca l high- 
latitude dust observed by DIRBE pwek et alj 119971 : lArendt et all [1998;). 
Right: Linear polarization vs. A for emission perpendicular to the local mag- 
netic field. 



are able to find size distributions and alignment functions (degree of 
alignment as a function of grain size) that are consistent with observations of in- 
terstellar extinction and starlight polarization. Figure IHI/e/t) shows the infrared 
emission from their models when the dust is heated by the local interstellar radi- 
ation field. The models are all in good agreement with the far-inf rared emission 



per H observed by DIRBE (|Dwek et al.lll997l : lArendt et al.lll998l ). 



Fig. \^right) shows that the different models predict quite different degrees 
of polarization, and with different dependences on wavelength. This is good news 
- it means that upcoming multiwavelength measurements of linear polarization 
by Planck should be able to reject some (or perhaps all!) of these models. 



5. Where Do Interstellar Grains Come Prom? 

The evolution of interstellar dust involves a complex interplay of many poorly- 
understood processes, in an arena - the ISM - that we also have a very incom- 
plete picture of. Ab initio approaches to the evolution of interstellar dust are 
therefore premature. 

Under these circumstances, the best approach to understanding the evo- 
lution of interstellar dust is to first determine what kind of dust is present in 
the ISM today. We can then try to construct an evolutionary scenario that is 
consistent with the grains that are observed to be present today. 
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Figure 5. Baryon budget for the Milky Way disk (see also Table 1). 



To summarize: Observations of interstellar gas-phase abundances ("deple- 
tions"), extinction, polarization, and emission by interstellar dust, and models 
that are able to reproduce those observations, tell us that: 

1. The abundant elements Mg, Si, Fe reside primarily in dust. 

2. Perhaps 2/3 of C is in dust. 

3. The Si is predominantly in amorphous silicate material. Based on available 
abundances, the composition is likely to be (approximately) Mgi.iFeo.9Si04. 

4. There is a substantial population of PAHs that contains ~ 10-20% of the 
interstellar C in the Milky Way; ~4-5% of the total grain mass is con- 
tributed by PAHs i n the Milky Way a nd other star-forming galaxies with 
~solar metallicity ( Draine et al.l [20071 ). 

Whatever the processes of interstellar grain evolution, they must be able to 
account for the above facts. 



5.1. Mass Budget for the Milky Way ISM 

Figure [5] illustrates the overall flow of mass into and out of the ISM. The mass 
of interstellar gas in the Milky Way is ~ 5 x 10^ Mq, divided approximately 
equally between atomic and molecular gas. Gas is being converted into stars 
at a rate of ~ 3M0/yr. If there were no source of gas, the ISM would be 
consumed on a relatively short timescale ~ 5 x 10^Mq/(3 MQyr"^) ~ 1.7 x 
10^ yr. However, there is return of material to the ISM from stars via various 



10 



B. T. Draine 



kinds of stellar outflows (stellar winds, planetary nebulae, novae, supernovae) 
and, in addition, there is evidence for infall of low metallicity gas. The rates are 
uncertain, particularly for infall, but it appears (see Figure [5|) that the net rate 
at which the ISM is losing mass is ~ IMq/jt, so that the characteristic decay 
time is of order 5 x 10^ yr, which seems a plausible number for a galaxy that 
has b een making stars in the disk at a more-or-less steady rate over the last ~9 
Gyr (|Rocha-Pinto et al.ll200d ^. 



6. Production of Stardust 

Some of the stellar outflows returning matter to the ISM are dusty. Outflows 
from cool oxygen-rich red giants contain silicate grains, and the outflows from 
cool carbon stars contain solid carbonaceous grains, in some cases with SiC as 
one of the carbon-bearing materials. The ejecta that form planetary nebulae are 
generally dusty. Some novae produce dust, but the net contribution from novae 
appears to be small. 

At least some Stardust comes from supernovae: supernova-produced dust 
grains have been identified in meteorites, and the internal extinction and IR 
emission from some young supernova remnants is evide nce for dust formation . 
Supernovae have even been described as "dust factories" ( Sugerman et al.ll2006l ). 



However, the average mass of dust per SN do es not appear to be la rge: the Type 
II SN 1987a produced ^ 8 x lO'^ Mq of dust (|Ercolano et al.ll2007l ^. The current 



record-holder for d ust production is the Type II SN 2003gd, which produced 
^ 0.02Mq of dust ( Sugerman et al.|[20 06). There is as yet no evidence of dust 



formation in Type la SNe. While only a very small number of SNe have as 
yet been studied, it appears that SNe produce, on average, ^ O.OIMq of dust, 
containing ^ 10% of the condensible elements. It should also be kept in mind 
that the dust that has been detected in these supernovae is moving at high 
speeds (~ 10^ km s~^), and is expected to undergo erosion by sputtering and 
possible destruction in grain-grain collisions when it is eventually decelerated to 
the low velocities of the quiescent ISM. While dust production unquestionably 
occurs, "dust factory" may be a misnomer. 



Table 1. Injection of Gas and Stardust from Stellar Sources 



gas 


dust 


Stellar Source 


0.4 
0.5 
0.06 


0.002 
0.0025 
<0.0001? 


Planetary Nebulae (— 0.3/yr) 
Red Giant, Red Supergiant, C star winds 
OB, WR, other warm/hot star winds 
SNe (1/100 yr, ~ IO'^Mq dust/SN?) 
Novae (100/yr, IO-'^Mq dust/nova?) 


0.25 
0.01 


0.0002? 
0.00001 


-1.2 


-0.005 


All stellar sources 



An accurate inventory of the different sources of Stardust is not available, 
but provisional estimates are offered in Table 1. The overall gas return rates 
from the different stellar sources in Table 1 can be estimated reasonably reliably. 
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given our knowledge of the IMF ( Kroupal l200ll ): (1) stars with initial masses 



above ~ 8Mq become supernovae, leaving behind either neutron stars or black 
holes. (2) Initial masses between O.QMq and ~ 8Mq leave behind white 
dwarfs of varying masses. (3) Stars with initial masses below ~ 0.9Mq do not 
evolve on the timescales of interest. 

We do not have a good understanding of either the stellar wind phenomenon 
itself, or the formation of dust in some stellar winds, so the estimates of dust 
injection are less reliable. The estimates for planetary nebulae and winds from 
cool stars assume that the dust mass is 0.5% of the gas mass, which corresponds 
to efficient condensation of the condensible elements. According to the estimates 
in Table 1, Stardust is entering the ISM at a rate Mstardust ~ O.OO5M0/yr. 

7. Stardust Lifetimes 

Newly-made Stardust moves away from its stellar source and begins a journey 
through the ISM. If nothing else happens, the grain will eventually find itself 
in a parcel of gas collapsing to form a new star. If this were the only process 
that removed grains from the ISM, the grain lifetime would be the lifetime of 
material against star formation: ra = rgp ~ Mism/SFR 2 x 10^ yr. 

The ISM is a violent neighborhood, agitated by fast stellar winds and super- 
nova explosions. A grain that is overrun by a fast shock is subject to sputtering 
by the sho cked gas. Sputtering in sh ock waves has been studied by a number 
of authors. iDraine Sz Salpeter ( 19791 ) found that a = 0.1/xm silicate grains had 



> 50% of the grain material returned to the gas in radiative shocks with shock 
speeds Vg > 200kms~^. 

We can estimate the global effects of SNe by a very simple argument. A 
supernova blastwave with initial energy £"0 = lO^^erg in gas with density rin ~ 
lcm~^ is Sedov-like until the shock speed drops to ~ 200kms~^. During the 
Sedov phase, 

Mv^g^Eo , (2) 
and therefore the mass of gas shocked at shock speeds > 200kms~^ is 

M(.. > 200kms-i) « ^ « 126OM0 . (3) 

Thus one supernova can "clean" ~ 1000 Mq of ISM of dust. If Eq = lO^^erg is a 
representative value for the kinetic energy of a SN, and the SN rate in the Milky 
Way disk is ~ 1/(100 yr), then the probability per unit time of an interstellar 
grain being overtaken hy a Vg > 200kms~^ shock is 

„i _ 1260MQ/100yr _ 1 

^ 5 X 109 M0 4 X lOSyr ' ^ ' 

There are other destructive processes as well, particularly grain-grain colli- 
sions in slower shock waves . Detailed studies ( Barlow &: Silklll977l : lDraine Sz Salpeter 



I1979I : iDwek Scalolll980l : Ijones et al.lll994l ^ that try to take into account the 
complications of a multiphase ISM find solid matter lifetimes Td ~ 4 x lO^yr, 
consistent with the simple estimate above. 
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To restate this: for a randomly-selected Si atom in a grain, the probability 
per unit time that it will be returned to the gas is t^^ ~ 1/(4 x 10® yr). In 
steady state, Stardust injection = Stardust removal, and we can estimate the 
mass of Stardust that should be present in the ISM. 



-^Stardust ~ -^Stardust 



^SF + ^d ' 



(5) 



from which we estimate Mgtardust ~ 1-6 x lO^AfQ. However, the total mass of 
refractory interstellar dust is A/jsmdust ~ 0.005 x Mjsm ~ 2.5 x 10''Afsun- Thus 
we conclude that 

Mstardust _ 1.6 X 10^ Mq _ _ 

Afismdust ^ 2.5 X WMq - ■ ■ U 



Most interstellar dust is not Stardust. 



Stardust accounts for only ~ 4% of the total mass of interstellar dust. 

Another way to approach the problem is to estimate the gas-phase abun- 
dance of Si. Let /^o(Si) be the fraction of Si entering the ISM from stellar 
sources (stellar winds, planetary nebulae, super novae...) that is in solid form. 
Because supernovae a ppear to condense ^ 10% of their condensibles into grains 
( Ercolano et al.l [20071 ) and the winds from many stars (e.g., O stars) appear to 



be dustless, we may estimate /*o(Si) < 0.5. The fraction of interstellar Si that 
is in Stardust will be 

/'^(Si) = Y^^^«0.2Ao(Si)«0.1 . (7) 

But in fact interstellar Si is almost always found to be strongly depleted, with 
/(Si) ~ 0.9 in dust. Therefore, we reach the same conclusion as before: most of 
the Si in interstellar grains is not in Stardust - it must have been added to the 
grains in the ISM. 

This is not a new result - iDraine &: Salpeter ( 19791 ) reached the conclusion 



that grain destruction was rapid and that regrowth of dust in the ISM was 
required to explain the observed depletions. The numbers basically haven't 
changed appre ciably since then; the argument has been reiterated a number 
of times (e.g. , iDwek fc Scalol Il980l : iDraind Il990l : IWeingartner fc Draind Il999l : 



Draind[2006bl ). Nevertheless, some authors continued to hold the view that the 



solids in the ISM were primarily formed in stars. 

For example, a number of papers have argued that the low abundance of 
crystalline silicate material in the ISM requires that crystalli ne Stardust materia l 
must be amorphized by cosmic ray damage in the ISM (e.g. iBringa et al.|[2007l ). 
However, the low abundance of crystalline silicates in the ISM can be easily 
understood: (1) grain destruction reduces the abundance of Stardust in the 
ISM, and (2) the silicate material grown in the ISM is amorphous. 

Isotopically anomalous silicates found in interplanetary dust particles (IDPs) 
are presumably Stardust. Techniqu es have recently been developed for study of 
Stardust silicates in meteorites (e.g. iNguven fc Zinneij[2003 : iKeller fc Messenger! 
20081 ). Recent studies of Stardust silicates identified in this way find that it is 
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70 amorphous, and ~20% crystalline material (L. P. Keller 2008, talk given 
at this conference). If this 4:1 amorphousxrystalline ratio applies to the overall 
production of silicate Stardust, then the observed u pper limit of ~2% on the 
crystalline fraction in the ISM ( Kemper et al.l [20051 ) implies an upper limit of 
2% X (4 + 1) = 10% on the fraction of interstellar silicate material contributed 
by silicate Stardust. This is in agreement with our estimate in eq. ([7]). The 
low crystalline fraction in interstellar silicates does not require that some of the 
crystalline material be amorphized - the expected destruction processes in the 
ISM suffice to keep the crystalline fraction below the current upper limit. 



8. Growth of Solid Material in the ISM 

If most interstellar solid matter is not Stardust, it follows that 



Most of the material in interstellar grains was formed in the ISM. 



The growth of grain material takes place on the surfaces of already-present grains 
- a process also referred to as "depletion" of condensible elements from the gas. 
We know that this process takes place, because we observe sightline-to-sightline 
variations in the gas phase abundances of the "depletable" speciesQ 
The depletion process has three stages: 

1. An atom or ion of the element must collide with a preexisting grain surface. 

2. The atom must become bound to the grain material in a way that allows it 
to be retained in the presence of fluxes of ultraviolet radiation and reactive 
species such as atomic H and O. 

3. The resulting grain material(s) will undergo heavy UV irradiation - what- 
ever survives (or is produced by) UV photolysis will be what is present 
in the ISM today. Cosmic ray irradiation also occurs, but is of secondary 
importance. 

Step 1 is just a matter of kinetics: The rate for a given ion to hit a grain will 
be determined by the total geometric area provided by the grains - the surface 
area will be dominated by the smallest grains - and by the random velocities of 
atoms and ions and the (possibly large) velocities of grains relative to the gafl 
In diffuse regions (where UV radiation is present) the ion-grain collision rate 
will be affected strongly by the charge state of the grains. Collisions between 
positively charged grains and ions such as Si"*" will be suppressed, but part of 
the grain population will be neutral, and the smaller PAHs are expected to have 



^ Mg, Al, Si, P, CI, Ca, Ti, Cr, Fe, and Ni are the prime examples of elements whose gas-phase 
abundances show order-of-magnitude variations from one sightline to another, the result of 
depletion onto grains. 

^ lYan et al] l|2004l 'l conclude that MHD turbulence will cause a > 0.04/im silicate grains in the 
"cold neutral medium" (tih ~ 30cm~^, T ^ 100 K) to move relative to the gas with velocities 
^ Ikm/s, large compared to the mean thermal speed {SkT /■Km)^^'^ — 0.27km/s of a Si ion. 
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an appreciable fraction that are negative, with which ions will have enhanced 
collision rates. 

Let Sd,2i ■ lO^^^cm^/H = dust geometric cross section per H nucleon. From 
modeling the UV extinction, we know that Sd,2i ^ 1- Let Va be the velocity of 
an atom relative to the grains. The timescale t^cc for this atom to collide with 
a grain surface is 

race = — = 1.0 X lO^r . (8) 

This time is relatively short. In the case of ions, the accretion rate will be 
modifie d by Coulomb repulsion (b y positive grains) or attraction (by negative 
grains). IWeingartner &: Draind (jl999) discussed the kinetics of accretion of ions 
(Ti+ was used as an example) onto grains in diffuse interstellar clouds, with 
attention to the charge distribution of the PAHs. They found that the very 
large surface area of PAHs resulted in a depletion rate 

race ~ 2 X loVr . (9) 

The fact that this is a factor 50 smaller than eq. (8) can be attributed to two 
factors. First, the total surface area in small particles is considerably larger than 
10~^^cm^/H. Second, carbonaceous grains smaller than 10~^ cm and silicate 
grains smaller than 2 x 10~^ cm were negatively charged, with Coulomb focussing 
increasing the collision rates by large factors. Evidently depletion and grain 
growth can take place on time scales ~ 2 x lO^yr in ordinary HI clouds, with 
nu ~ 30cm~^. 

We conclude that the low gas-phase abundances of elements such as Si and 
Ti can be understood based on rapid rates of depletion from the gas phase for 
the conditions found in ordinary HI clouds - the "cold neutral medium". The 
bulk of interstellar grain material is the result of accretion from the gas in the 
cold ISM. 



9. What Materials Form? Why Do Some Elements Not Deplete? 
9.1. Thermal Desorption 

In typical diffuse clouds, some elements (e.g., Ca, Al, Si, Fe) are heavily "de- 
pleted" fror n the ga.s, whi le others (e.g., S, Cu, Zn) show much less of a tendency 
to deplete ( Jenkind 20041 ) . To try to understand how this selectivity can occur. 



let's examine the steps in the depletion process. 

Grains with radii a ^ O.Ol^m remain very cold in the diffuse ISM, with 
temperatures T ^ 20K. At these temperatures, almost any atom or ion imping- 
ing with a kinetic energy ~ O.OleV is expected to have a high probability for 
"sticking" to the surface at least temporarily (where "temporarily" means long 
compared to the characteristic vibrational period ~ 10~^^s), with the exception 
of the inert gases such as He and Ne. Normal physisorption via van der Waals 
forces will bind species such as S, Cu, and Zn to grain surfaces with binding ener- 
gies B ~ 0.05 - 0.2eV. Thermal desorption will take place with a probability per 
unit time ~ vq ex.p{— B / kT) , with ~ lO^^s"^. For B = 0.05 eV and a grain 
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temperature T = 20K, the lifetime against thermal desorption would be ~ 1 s, 
whereas if B = 0.10 eV, the lifetime becomes 5 x 10^ yr. Therefore if an atom 
is to remain on the grain surface long enough to matter, it must have a binding 
energy B > 0.1 cV. Unfortunately, the binding energies of the atoms of interest 
(C, Mg, Si, Fe) to surfaces of interest (amorphous silicate or carbonaceous ma- 
terials) are not known, but it is plausible to imagine that those elements that 
do deplete are held to the grains surface by binding energies B > 0.1 eV. 

9.2. Hydrogenation 

We must now estimate some other rates. Consider one "surface site" on a 
grain, where some species X is adsorbed. The area of this surface site will be 
A « 10~^^cm^. This surface site is being bombarded by H atoms at a rate 

,,,, / kT \^/^ ^ If n(H) \ / r 

The impinging H atoms might form chemical bonds with the adsorbate X; 
whether they will do so is not known (for some X there may be an energy 
barrier in the "reaction coordinate" necessary to produce XH, preventing re- 
action at low energies) but it does seem possible that incident H might result 
in hydrogenation of some physisorbed atoms, such as S. The resulting molecule 
might be ejected from the grain surface, or might remain there. 

In H I clouds, O atoms are the third most abundant species in the gas (after 
H and He), and will collide with a surface site at a rate 

"^"^^2^) ^ " 3 X IQSyr (sOcm-^) (iOOk) ' ^^^^ 

9.3. Photodesorption or Photolysis by UV Radiation 

Dust grains in diffuse clouds are exposed to diffuse starlight, including ultraviolet 
radiation from hot stars. The starlight intensity is such that, for example, a C 
atom in the gas (with ionization energy 11.26 eV) will be photoionizcd at a 
rate <^(C C+) = 3 X 10""^'' s~"^. Species on grain surfaces will in general have 
allowed transitions to excited states with E < 13.6 eV, and these transitions will 
also be excited by the starlight background, and we may assume a characteristic 
rate for UV excitation 

C« 10-i°s-i = — ^ . (12) 
^ 300yr ^ ' 

What happens when the adsorbed species (which might be an atom X, a hydride 

XH, an oxide XO, a hydroxide XOH, ...) absorbs a UV photon and makes a 
transition to an electronically excited state? The answer is: it depends. Some 
materials (e.g., stainless steel or quartz) can undergo heavy ultraviolet irradia- 
tion without being affected. Other materials (e.g., polyethylene) are chemically 
altered by UV irradiation. Some physisorbed or chemisorbed species will un- 
dergo photoexcitation to an electronic excited state that may actually be repul- 
sive with respect to the surface, resulting in ejection in a time ~ 10~^^ s after 
being photoexcited. Laboratory studies of H2O and mixed ices irradiated by 
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UV radiation ( Westlev et al.|[l995l : lOberg et al.ll2009al lbl) show that the surface 
erodes at a rate consistent with ejection of molecules from the surface monolayer 
with a high yield per photoexcitation in the surface monolayer. 

Comparison of eq. (10-12) shows that a species X residing on a grain surface 
will have a good chance of encountering an impinging H atom before anything 
interesting happens, but will be very unlikely to encounter any other impinging 
species (e.g, C, O, Si) before undergoing photoexcitation. The adsorbed species 
will undergo ~ 10^ photoexcitations before a monolayer of material (other than 
adsorbed H) can be deposited on top of it. 

Because photoexcitation happens so frequently, the only materials that can 
form by accretion in the diffuse ISM are those that are not destroyed or altered by 
ultraviolet irradiation - this is what selects for amorphous silicate material and 
carbonaceous material as the two condensates. There may be other materials 
present, of course - e.g., Fe oxides, or even metallic Fe - but observations of 
interstellar dust appear to be consistent with amorphous silicates containing 
most of the condensed material that is not carbonaceous. 

Ultraviolet irradiation may also explain how it is possible for the ISM to be 
able to grow two distinct grain materials - carbonaceous material and amorphous 
silicate material - out of a single gas mixture. One can imagine the following 
scenario: Suppose the grain population already has some amorphous silicate 
material exposed to the interstellar gas. When Mg, Si, Fe, and O atoms and 
ions arrive at the grain surface, they are able to grow additional amorphous 
silicate, adding one atom at a time. It may actually help to have ultraviolet 
radiation present, as the electronic excitations may allow chemical bonds to be 
rearranged to form amorphous silicate material as new O, Mg, Si, Fe atoms are 
added to the surface layer. 

On the other hand, what happens when a C atom, for example, lands on 
the amorphous silicate surface? One can imagine that the C atom physisorbed 
on the amorphous silicate surface might undergo photoexcitation to an excited 
state that is repulsive, ejecting it from the surface. Or perhaps the C would 
become hydrogenated or oxidized, with the resulting CH or CO undergoing 
photodesorption from the surface. Such processes could keep the amorphous 
silicate carbon-free. 

Similar processes may occur on exposed carbonaceous surfaces: impinging 
C atoms could grow new carbonaceous material, whereas impinging Mg, Si, Fe, 
etc. could be removed by some combination of reaction with impinging H or O, 
and photoexcitation by UV. 

In this scenario, we can imagine the carbonaceous material and amorphous 
silicate material being in separate grains, or perhaps even as separate regions 
of a single grain. In the latter case, one would expect that when grain-grain 
collisions cause grains to shatter, the fragmentation would preferentially occur 
along interfaces between the two materials, leading to "pure" fragments. 

Perhaps, then, the interstellar grain population actually is dominated by two 
distinct grain types: amorphous silicate grains and some form of carbonaceous 
material. 
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10. Growth of Dust at High z: the example of J114816+525150 

Observations of quasars and luminous gal axies at high reds hift have detected 
large masses of dust in a number of systems ( Wang et al.ll2008l ). The prime exam- 



ple is the galaxy associated with the z = 6.42 QSO SDSS J114816+525150 (here 
after J1148+5251), with estimated dust mass Mdust ranging from 2 x IO^Mq 
(jPwek et alJl2007l i to 7 x W^Mq (iBertoldi et alJl20C)3al ^. Molecula r gas appears 



to account for a large fraction of the ~ 5 x W^^Mq dynamical mass (jWalter et al 
l2004l ^. With Mgas < 5 x 10^^ Mq, and Mdust 2 x W^Mq, this system has a 
dust/gas mass ratio Mdust /-Mgas > 0.004 - comparable to or exceeding the Milky 
Way value ~0.005. 

The time available for stellar evolution prior to z = 6.42 is limited: if star 
formation began at z = 10, the oldest stars are only 400 Myr old at z = 6.42, and 
only massive stars will have been able to evolve ~ low-mass stars have insufficient 
time to evolve to the asymptotic giant branch which dominates production of 
Stardust in the Milky Way. This has led a number of authors t o propose that 



supe r novae are responsible for the dust in high- 2: gala x ies fe.g. iMaiolino et al 
20041 : ISugerman et al.l[200i : lBianchi SchneideJl2007l l. iDwek et al.l (|2007l 1 dis- 



cuss the dust in J1148+5251 and conclude that supernovae would have to pro- 
duce ^ IMq of dust per supernova to explain the observations, but note that 
this is considerably in excess of what has been observed in SN ejecta. 

J1148+5251 con tains a large mass of rnolecular gas, detect ed in CO J = 7 — 
6, 6 ^ 5, and 3 ^ 2 (jBertoldi et al.ll2003bl : IWalter et al1l2004l '). The CO that is 



observed, and the H2 that must accompany the CO, is not supernova-produced: 
even if those molecules do form in the ejecta, they are efficiently destroyed by 
the reverse shock when the high-velocity ejecta are decelerated. 

Instead, the H2 must be formed by the same mechanism that dominates 
H2 formation in the Milky Way: catalysis on grain surfaces. Given that the 
massive stars present in these galaxies will destroy H2 molecules - primarily 
through photodissociation - each H nucleon in the gas must, on average, have 
collided with grain surfaces many times in the age of this galaxy. Metal atoms 
and ions move only a few times more slowly than H, and will also collide with 
grain surfaces many times; if they stick, they will form new grain material. As 
discussed above, this is the process that dominates grain formation in the Milky 
Way, and there is no reason not to expect it to formation of grain material in 
J1148+5251. 

Supernovae are of course required to produce the metals that compose the 
grains, and to provide some supernova-condensed "stardust" to provide "seed" 
surface area on which to grow more material in the ISM, but the bulk of the dust 
mass in high-z galaxies appears likely to be the result of grain growth competing 
successfully with grain destruction in the ISM. 



11. Summary 

The principal points of this paper are as follows: 



1. The principal observational constraints on models for the interstellar grain 
population are summarized. The observed wavelength dependence of ex- 
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tinction (Fig. 1) and polarization of starlight, together with elemental 
abundances (total and observed in the gas) continue to provide the strongest 
constraints on what interstellar dust can be. Additional information is pro- 
vided by infrared emission features at 3.3, 6.2, 7.7, 8.6, 11.3, 12.0/xm that 
are identified as coming from PAHs following single-photon heating. 

2. A grain model consisting of a population of amorphous silicate spheres 
and a population of carbonaceous grains, including PAHs, can reproduce 
the observed wavelength-dependent extinction by dust. This model also 
successfully reproduces the observed infrared emission from interstellar 
dust (Figs. 2, 3). 

3. If the amorphous silicate grains and carbonaceous grains are given spheroidal 
shapes, the observed polarization of starlight can also be reproduced. 

4. The infrared emission spectra calculated for these spheroidal grain mod- 
els (Fig. 4) are consistent with the observed infrared emission from high- 
latitude regions in the Milky Way. The models predict different degrees of 
polarization as a function of wavelength - observations with Planck will 
test these predictions. 

5. The evolution of interstellar grains in the Milky Way is discussed. Grain 
destruction in the ISM is such that ;^10% of the interstellar dust mass 
consists of "stardust" from stellar sources, including supernovae. The bulk 
of interstellar dust has been grown in the ISM. 

6. The amorphous silicate and carbonaceous materials present in interstellar 
grains are the product of grain growth in the presence of strong ultraviolet 
radiation, that will photoexcite and possibly photodesorb species from 
the grain surface. The ultraviolet radiation may be key to maintaining 
separate populations of silicate and carbonaceous grains even though most 
grain growth occurs out of an average-composition ISM. 

7. Dust in high-redshift systems such as Jl 148+5251 must be predominantly 
the result of grain growth in interstellar clouds. A small amount of dust 
from supernovae is sufficient to initiate grain growth in the ISM. 
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